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Abstract 

Let us call "Gibbs energy" the quadratic form occurring in the maxi- 
mum likelihood (ML) criterion when fitting a zero-mean multidimensional 
Gaussian distribution to one realization. We consider a continuous-time 
Gaussian process Z which belongs to the Matern family with known regu- 
larity index v > 1/2. For estimating the range and the variance of Z from 
noisy observations on a dense regular grid, we propose two simple estimat- 
ing functions based on the conditional Gibbs energy mean (CGEM) and the 
empirical variance (EV). We show that the ratio of the large sample mean 
square error of the CGEM-EV estimate of the range-parameter over the one 
of its ML estimate, and the analog ratio for the variance-parameter, both 
converge (when the grid-step tends to 0) toward a constant, only function of 
z/, surprisingly close to 1 provided v is not too large. Extensions of this ap- 
proach, which may enjoy a very easy numerical implementation, are briefly 
discussed. 

1. Introduction 

We consider time-series of length n obtained by observing on a dense regular 
grid, a continuous-time process Z which is Gaussian, has mean zero and an au- 
tocorrelation function which belongs to the Matern family with regularity index 
v > 1/2. This family is commonly used, e.g. in geostatistics (see Stein 1999); recall 
that v = 1/2 correspond to the well known exponential autocorrelation (in other 
words, Z is a stationary Ornstein-Uhlenbeck (O.U.) process). The general defini- 
tion of Matern processes can be easily formulated in terms of the Fourier transform 
of their autocorrelation function, namely the spectral density over (— oo, +oo): 
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ft,bfii u ) = b 9t,e( u )' where 9l,ei u ) = (m - rrr~~ - 

In this paper the constant c u = -j=^l is chosen so that J^g^^^du = 1. Thus 
6 is the variance of Z(t) and 6? is the so-called "inverse-range parameter" (in fact, 
it is Iv 1 ! 2 /6 which can be interpreted as the effective range or "correlation length", 
cf Stein (1999, Section 2.10); we will often drop the term "inverse". 

We are concerned here with dense grid in the sense that the interval 5 between 
two succesive observations is small relatively to 1/9. The considered processes 
being mean square continuous, this means that two successive ideal (i.e. with- 
out measurement error) observations are "strongly" correlated. See Stein (1999, 
Chapter 3), for such a setting which shows that standard large-n asymptotic anal- 
ysis followed by a (less standard) small-5 analysis yields useful insights and good 
approximations for various real (finite-size) problems. 

In this article, we assume that v is known and we mainly study settings where 
there are Gaussian i.i.d. measurement errors, or, equivalently for the paramet- 
ric inference point of view we take here, there is a so-called nugget effect. No- 
tice that introducing such a nugget effect is actually not too restrictive since, of 
course, data always have a finite number of digits, so one may easily accept to 
include, in the parametric model, a known nugget equal to the squared finite pre- 
cision; furthermore this is known to often remedy to numerical instabilities due 
to ill-conditionning of the linear systems we may encounter in the no-nugget case. 
However, in this first study, we restrict ourselves to the case where the suspected 
measurement errors have a known variance. We thus assume that (possibly after 
appropriate rescaling) the observations is a vector y of size n whose conditional 
law given Z is 

y|z ~ N(z, I) where z ~ N(0, boRe ), 

with / denoting the identity matrix and Re the autocorrelation matrix of the col- 
umn vector z = (Z(S), Z(25), . . . , Z(nS)) T , i.e. the Toeplitz matrix with coefficients 
given by 

/oo 
-oo 

Recall that expressions in terms of Bessel functions are well known for K v> q (see e.g. 
Rasmussen and Williams (2006) for explicit very simple expressions for v = 3/2 
and 5/2). 

It is often of great interest to be able to "effectively reduce" the number of 
parameters, especially when computing the likelihood function is costly. This is 
classically done in the no-nugget case (i.e. bRe + 1 is replaced by bRe, and y by z, 
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in the likelihood function) by noticing that the maximizer b of the likelihood for 
any fixed 9 is simply: 

b(0) = (l/n)z T R e - 1 z. (1.2) 

See Zhang (2004) also for numerical improvements produced by such a "concen- 
tration" of the likelihood. Zhang and Zimmerman (2007) recently proposed to use 
the classical weigthed least square method (not statistically efficient but much less 
costly than maximum likelihood (ML)) to estimate the range parameters, next, to 
plug-in these parameters (the 9 here) in the likelihood which is then maximized 
with respect to b (the solution, say 6ml {9), being the explicit (1.2) in the no-nugget 
case, obtained iteratively by Fisher scoring otherwise). The idea underlying this 
method is that, at least for the infill asymptotics context, even if 9 is fixed at a 
wrong value 6*i, the product 6ml {9i)9\ v still remains an efficient estimator of b 9l u 
(see Du et al. (2009) for recent results of this type). 

The method we propose here, firstly reverses the roles of variance and range, 
in that it is based on a very simple estimate for the variance, namely the empirical 
variance in the no-nugget case, and its corrected version for biais otherwise, which 
is simply defined by 6 EV := n~ 1 y T y — 1. Secondly we propose to replace the 
maximization of the likelihood w.r.t. 9 by the very simple following estimating 
equation in 9, in the with-nugget case: solve, with b fixed at &ev 

y T A bfi (I - A b:9 ) y = trA bfi where A bfi = bR e (I + bRo)' 1 . (1.3) 

In the no-nugget case, this equation in 9 is simply replaced by z T Rg^ 1 z = n&EV- 
One may call "Gibbs energy" (GE) of the underlying process the quantity z T R e ~ 1 z 
and it is easily seen that b (y T A b ^ e (I — A bt0 ) y + tr(J — A btd j) is the conditional 
Gibbs energy mean (CGEM) obtained by taking the expectation of z t Rq~ 1 z, con- 
ditional on y, for the candidate parameters b, 9. So equation (1.3) in 9 will be 
called the CGEM-EV estimating equation (GE-EV in the no-nugget case) and we 
will denote by #geev this new range parameter estimate. 

Notice that it is easily checked that equation (1.3) is equivalent to equate to 
the derivative w.r.t. b (and not w.r.t. 9\) of the likelihood. This proposal is based 
on the following two ideas. First, since it is quite plausible that the above idea, 
underlying Zang and Zimmerman (2007) 's proposal, remains true for a random 
9, then, instead of "fixing" 9 at #1, one may as well adjust 9 so that &ml(#) 
coincides with a given value bi for the variance; and, denoting 9\ the so-obtained 
9 (thus (b±,9i) must cancel the 6-derivative of the likelihood), the product bi9\ u 
will plausibly be an efficient estimator of &o$o^- Second, it is known that, at 
least in the no-nugget case, the moment Z(5j) 2 yields a successfull estimating- 
equation in the case v = 1/2 to estimate 9 when b 9Q U * s known (see Kessler 
(2000)). So it is tempting to select bi := b E y, the emprical variance estimate of 
b . The second point is an admittedly weak justification of this method. However, 
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notice that Yadrenko (Chapter 4 Section 3) has studied, for quite general, mean 
square continous, stationnary isotropic random fields the properties of a continuous 
version of the empirical variance in the periodic case, i.e. the Lebesgue integral 
over sphere of Z(t) 2 , and he has listed appealing properties of such estimates of 
the variance of the field. The theorical results we give in this article, will provide 
in our context, a much stronger, and rather unexpected, justification, for v not too 
large (which is the case in numerous applications, see Stein (1999), Rasmussen et 
al. (2006)). 

Of course, in our time-series setting, there now exist rather good implemen- 
tations of ML; see Chen et al. (2006) for the case when the correlation is rather 
strong. However, our objective here is to provide first insights into the capability 
of this method for more computationally complex settings. Indeed, this approach 
is not limited to observations on one dimensional lattice, and is potentialy not lim- 
ited to regular grid (a weighted version, with Rieman-sum type coefficients, of the 
empirical variance shoud then be used instead) and to homogenous measurement 
errors: some successfull applications for two-dimensional Matern random fields, 
using a randomized-trace approximation to trA^, are described in Girard (2009), 
along with some theoretical properties. This approach might be, in principle, 
applied to other two-parameters models than the Matern family (of course it is 
presently restricted to scalar 6). Encouraging experimental results are obtained in 
Girard (2009) for the commonly used spherical autocovariance functions. 

The rest of this article is structured as follows. Section 2 lists additional no- 
tations. In the asymptotic framework we adopt here (which may be thought as 
intermediate between the infill and increasing-domain frameworks) we first show, 
in Section 3, that, even when b is arbitrarily fixed to b 1: the CGEM estimating 
equation is a quite well-behaved estimating equation: it converges toward a mono- 
tonic equation whose the root Q\ satisfies b\Q\ v = b o 0Q U . Next we study in Section 4 
what is sacrificed when plugging-in the simple &ev m the CGEM estimating equa- 
tion, as compared to ML. The performance- loss is classically quantified by the 
asymptotic mean square errors of the two components of the CGEM-EV estimate 
as related to those of the ML estimate. We show that what is sacrificed, which 
depends on z/, is quite reasonable provided v is not too large, independently of bo 
and #o- Indeed asymptotic efficiency is reached as v decreases to 1/2. Proofs are 
outlined Section 5. 

2. Further notations 

Of course "time" could be replaced everywhere by "space": we choose this 
vocabulary since we use in several places of the paper the classical time-series 
theory. Thus we assume everywhere without loss of generality that defined by 
Z$(i) = Z(5i) is observed at times i — 1, 2, • • • , n. The spectral density on (—it, 7r] 
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of Z 5 is thus f vhfi = bg s ve with 



<*(•) : = s E ( lJ r Z ) = £ C (• + 2A;7r ) where a = se > 

k=—oo ' k=—oo 

the equality between the two sums resulting from the particular "variance-scale" 
parameterisation of the Matern family. We will also denote simply by g Ut0l (•) 
the second sum (a is the range-parameter for the series on Z). Simple closed 
expressions for g & u6 are available only when v — 1/2 is a small integer, says 0, 1 or 
2 : they then coincide with particular constrained ARMA spectral densities. 

Let us define the "Bayesian filter" a & be (-) and its "un- aliased" version a\ {•), 
for a = 59, by 



6 f \ -Jl/,0 \ / * I \ 

a b,ev> -~ uZS ( \ , To \-1 ' a i>,aU • 



^(•) + (2vr)- 15 bgtj-) + (2n)~^ 

As is well known af e is the spectral characteristic of the "n = +00" version of 
A b e (a convolution of series defined on Z) so that, e.g., lim n _ ++00 n _1 trA 6 = 
(2ff)- 1 /" w a^(A)dA for any fixed 5, b, 9. 

A function which will play an important role in this article is the derivative of 
\og{g 5 ve {-)) w.r.t. 9, and its un-aliased approximation (up to a factor 5): 

h% := d]ag(gl t e)/de, := d\og{g^ a )/da. 

For any / : [— tt, tt] — > R, s.t. /^[(^(A^/MdA 7^ 0, we define the weighted 
coefficient of variation of / by 

T- N If- f w f)Y -r-fwf 2 

T ( f\ J W L \J wJ \ S w J 1 U / J \2 

A&M/) : = 7 ^2 — = t \2 - x > where w := (°W ■ 



Above and throughout this paper, "J" will denote integrals over [— tt, tt]. We will 
also use the notation g (resp. h ) for the function g 5 ug (resp. h 5 ug ) when 9 = 9 
and also a := aj,^, and J (-) := Js,bo,e (-)- 



3. Some asymptotic properties of the CGEM estimating equation 

In this section, we established some consistency and identifiability properties 
enjoyed by the CGEM estimating equation in 9. They are merely encouraging pre- 
liminary results. We expect to give more complete results elsewhere. To simplify 
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our statement, we consider, the normalized version y T Ab i g(I — A^s)y /trA^g = 1 
which is, of course, an equivalent estimating-equation. 

For any positive 5, b, 9, bo, 0q, we can define the following weighted mean: 

Theorem 3.1. For any fixed b, 9, we have the following convergence in probability: 
y T A be (I — A bg )y 

'-— ■ 1 -> ip(5, b, 9, b , 6 ) as n -> oo, 

trA 6) e 

and ip has a very simple small-5 equivalent: 

(b 2u \ 



The first part of this theorem is in fact not restricted to the Matern fam- 
ily; this is a rather direct consequence of universal results of time-series the- 
ory about quadratic forms constructed from (product of power, possibly nega- 
tive, of) Toeplitz matrices. In order to be allowed to apply such classical results 
(for instance, those stated by Azencott and Dacunha-Castelle (1986) and used in 
their analysis of maximum likelihood estimation, or its Whittle approximation), 
it is sufficient that (9, A) i— > gl e {X) be three-times continuously differentiable on 
x [— 7r, 7r] where 6 is a compact interval not containing 0; and this can be checked 
by applying the classical Weierstrass M-test. 

The second part uses several appromixations to integrals which can be obtained 
by the technics of Section 5. 

Theoreme 2.1 thus says that the limit of the normalized CGEM equation has 
a small-5 equivalent which is a very simple monotonic function of b9 2u . Clearly if 
b is fixed at any value b\, then the unique root 9\ of this large-n-small-5 equivalent 
equation will satisfy b\9\ v = b 9 2u . This gives support to the first of the two 
underlying ideas given in the introduction. 



4. Mean square error inefficiencies of CGEM-EV to ML for the variance 
and range parameters 

Let (&mlj 9ml) be a maximizer of the likelihood function on B x 6 where B 
and are compact intervals not containing and such that (b , 9 ) is in the 
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interior of B x G. Then, since the classical ident inability and regularity conditions 
are well fulfilled, for any fixed 5 > 0, it is a well known result of times-series 
theory (Azencott and Dacunha-Castelle (1986)) that (&ml,#ml) is a.s. consistent 
and satisfy: 
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1/2 
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0~12 



0\ 2 cr 2 



where 



0~12 



(Joq/io) 2 ^o(/io) 



b 2 J alhl 
-b J a 2 Q ho 

Ul 



Note that the factor (/ alh Q )~ 2 J (/i ) _1 is the determinant (> since ^o( - ) 
cannot be a constant function) of the information matrix and has this form only if 
/ alh 7^ 0. And the same regularity conditions on our time-series Matern model 
fl b9 are also sufficient to show (by the usual Taylor series argument) the first part 
of the following theorem: 

Theorem 4.1. Let (&ev,#geev) be a consistent root of the CGEM-EV estimating 
equation (1.3). If f alh Q ^ then 
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1/2 



b E y 

#GEEV 



bo 
0o 
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V12 V 2 



where 



Vl = blfa 2 g$ 

vu = -b c (J alh ) , and c = J (g /al). 
v 2 = c J al(J alh ) 

we have f a\hoj f a 2 , ~~ * ^ u /^o> an d defining by Ig t 
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(resp. Isb e := v 2/°~ 2 ) the asymptotic mean square error inefficiency of CGEM- 
EV to ML for bo (resp. for 9 ), these 2 inefficiencies have the following common 
limit 



l 8,b ,e 



v^F /2z/ + i\ 2 r(z/+i) 2 r(2z/+i) 



2u 



r(z/) 2 r(2i/ + 1) 



as 



5^0, for ie{l,2}. 



The first rather surprising fact is that this large-n-small-5 inefficiency is not 
function of the underlying bo, nor of 9q. Secondly, Table 4.1 clearly demonstrates 
that the CGEM-EV estimates are asymptotically nearly efficient provided v is not 
too large. Asymptotic full-efficiency is reached for v close to 1/2. 
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1/2 


1 3/2 


2 5/2 
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7/2 
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1 


1.04093 § 


1.18596 jjg 


1.33174 


1716 

1225 


1.46727 



Table 4.1. Limit of asymptotic MSE inefficiencies P s bo 0q , i £ {1, 2}, for some "typical" values of ^ 



Of course it would be interesting to study how fine and large the grid must be 
in order that the CGEM-EV statistical performances be comparable to ML. One 
may guess that the signal-to-noise ratio must also be "not too small", even if its 
value does not matter asymptotically. Since the main improvements in numerical 
performance occur for multidimensional process (with randomized-traces used in- 
stead of the exact traces in (1.3)), a first Monte Carlo study has been done for 
two-dimensional random fields in Girard (2009). 

This approch may be extended in a simple way to the case of unknown vari- 
ance, says a 2 , of the measurement error (or unknown nugget effect). Indeed it is 
well known (Stein (1999, Section 6.2)) that the simple average of local squared 
differences like (y^ — Vi-i) 2 would yield here a consistent estimate a 2 , and this 
is easily generalized to multidimensional processes which are mean square contin- 
uous. Notice that the Bayesian filter is still A bi g except that b now denotes the 
signal-to-noise ratio s.t. var(Z) = bo 2 . Then the whole CGEM-EV approach can 
be applied after having replaced y by its "standardized" version y/a e . It would 
be thus very interesting to study whether similar near-efficiency results still hold 
or if a more elaborate estimate for o e is required. 

Concerning the case with known noise and signal variances, even if the as- 
sumption that bo were known is rather restrictive, it may be worth to point out 
the following very neat result: denoting by #ml t ne ML estimate with b fixed at 
b , its asymptotic behavior is classically obtained 

n 1/2 (§ml - 0o) ->t> N ^0,471^ a 2 h 2 

and denoting by §ge a consistent solution of the estimating-equation (1.3) with 
now b fixed at bo, we have 




Theorem 4.2. // / a 2 h ^ then 

n 1/2 (0geo " 0o) ->2> N ^0,4tt J a 2 Q a 2 h 

and the ratio of this asymptotic variance over the one of9 MLo , says Is,b ,e > satisfies: 

^6o, e „ = 1 + J o(M^l as 5^0. 
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Asymptotic (i.e. large-n-small-5) full-efficiency is thus now enjoyed by the 
CGEM estimating function for any value of v > 1/2. 

To finish, let us comment on the somewhat more natural &EV-based alterna- 
tive to ML which would be akin to the hybrid estimate proposed by Zhang and 
Zimmerman (2007) : choose #r- so as to maximize the likelihood function with &ev 
plugged-in. As to the computational aspects, #h m ay be much less easily com- 
puted than #geev (for example when an iterative solver is used for A^ey) since 
it requires a log-determinant instead of a trace. So it is worth to point out that 
again a neat result (which can be similarly obtained by the technics used to prove 
Theorem 4.1) holds true here: the large- n mean square error of #h over the one of 
#geev tends to 1 as 5 — > 0. 

5. Proofs 

In the following, c denotes a constant that may change from line to line (that 
is, c will only possibly depend on v and on the positive lower and upper bounds 
for the candidate 6's and those for 9 (says b,b,9_,9)). And 0(5 2u ) denotes a term 
whose absolute value is bounded by c5 2u . Without loss of generality we assume 
that S < 1. 

Firstly, it is easily checked from the definition of g£ ;Q ,(A), that the un-aliased 
version of h s v 9 has a very simple form: 

KJX)=a- 1 ^2u-(2u + l)G(^jy where G(u) = (5.1) 

We first collect intermediary results in a few lemmas which might be of interest 
for other studies of this small-5 Matern time-series model. 

Lemma 5.1. We have 

cb 2v < 9 ; t5e (\) ^ gl fi {\) ^ g^ e (X) + 0(5 2 n, 

and 

9:,so( X )KeW-SK M (\)\=0(5 2 n- 

Lemma 5.2. We have 

^(^i(A)aUA)+0(0, 
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and 



Proof of these Lemmas. From their definition, we have g v>a (X) ~ 9t, a ( X ) = 
J2k^o 9t,a( X + 2nk). For any A G [— 7r, 7r], for any k > 1, the monotonicity of g* a 
implies #* a (A + 2ttA;) < g^ a {2n(k - 1/2)) < ca 2u /(k - l/2) 2iy+1 ; and summing 
these terms over k = 1, 2, • • • , thus gives a term 0(a 2u ). This is shown similarly 
for the sum over k = — 1, — 2, • • • . We now prove the second part of Lemma 5.1. 
Omitting the index v and denoting g a := dg a /da, g* a := dg* a /da, one can decom- 
pose ag* a {5- l h 5 e - h* a ) = ag* a (g a /g a - g*/g* a ) = (g* - g a ) (ag a /g a ) + a(g a - g*). 
Observing that ag* a = 2vg* ua + cg* +la yields a(g a -g* a ) = 0(a 2u ). Lastly, 
we can bound a\g a /g a \ < a\g a /g* a \ < a\g* a /g* a \ + a\g a - g*\ / g* a < 0(1) + 
a \g a — g* a \ /ca 2u from equation (5.1) and the first inequality of Lemma 5.1. 

Lemma 5.3. If v > 1/2, then for any continuous function F : R — > R + such that 
< F(x)dx < +oo, we have, as 5 — > 

/ <oW F (jo) dX ~ 69 jZ F{x)(ix ' I \- a U x ^ F (je) dXr ^ 69 SZ F{x)(ix ' 

Proof. The first part of Lemma 5.2 and the boundedness of F permits us to 
replace af^ by the much more manageable a* b a (where a = 56) with an error 0(S 2u ) 
which is enough when v > 1/2 . Letting A = a e+ ( 2u /( 2u + 1 )) where e is any arbitrarily 
small constant > 0, we see from the monotonicity of g* , that inf A6 r_^i a£ a (A) > 

(l + (2ti /b) [^(A)] -1 )" 1 which tends to 1 since #* a (A) > a 2v / (2A 2 )^ +1/2 (from 

a small enough) which tends to +oo. And the same is true for [a^ a (A)] 2 . It 

remains to observe that both J^. F(^)dA and J_^F(^)dA are equivalent to the 
claimed term. 



Lemma 5.4. We have 

//»oo 
a^(A)dA ~ (50)^1 (2nc v b)^ J (1 + x^+^dx, 

J [al fi (\)] 2 d\~ (56)^i (27ic v b)^ J°° (l + x 2u+1 )- 2 dx. 

Proof. By Lemma 5.2 we can replace, with enough accuracy, the filter by its un- 
aliased version in these integrals. Next the claimed equivalents can be obtained by 
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2v 

a change of variable s = A/a 2 "+ 1 and an application of the dominated convergence 
theorem. 



We can now outline the proofs of Theorem 3.1 and Theorem 4.1, beginning by 
the second one. 



Proof of Theorem 4.1. i) First it is directly seen that i^/of = </o(^-o) c o- Let 
w(\) := [/^(A)] 2 . Consider first the constant c . Since l+J$ bt g(g^ d /w) is the ratio 
of fwf[gi/\)] 2 /w(\)d\, where [<£J 2 /w = [<f vfi + (2Kb)- 1 ] \ over [J^dA] 2 
( = 1), it suffices to observe that 

/ 9u, a 2 ~ J gtj < J \g»,* 2 - 9lJ\ < ca 2 » J (g„, a + £ >a ) < 2co> 

and that J g* v 2 ~ a^c 2 (1 + x 2s )~ 2v ~ x dx to see that this is also the dominant 
term of J8,b,e(gifi/ W )l f w. Secondly, the equivalent of Jo(/io) can be obtained as 
follows. The numerator of Js,bfi(h 5 v e ) is 



I 



w 



| w(h%-2u/9) 2 - 



I 



w 



w(hl fi - 2v/6) 



If we take for granted that the un-aliased Sh* se can be substituted for h 5 ve , then 
the numerator is simplified using equation (5.1), and its equivalent is deduced from 
Lemmas 5.3 (case v > 1/2) and 5.4 (with G a (X) := G(X/a) defined from (5.1)) as 
follows: 



' 1 w\(2u+l)GJ6\ 2 - 



w 



1 



w 



w(2v+l)G a /0 



((2i/+l)/0) 2 -^50 / G(x) 2 dx. 



I 



w 



Note that this can also be proved for v = 1/2 by using the closed expression of w 
available in this case. On the other hand, the denominator, also for the un-aliased 



which tends to (2v/6) 2 again 



version, is (J w) 1 / w(\)[2v - (2v + l)G(X/a)}/6 
by Lemmas 5.3 and 5.4. This proves the claimed result for the ^-inefficiency pro- 
vided the un-aliased substitution was authorized. To see this, e.g. for J w(h s u9 — 

2v/6) 2 , it suffices, to decompose (h 5 vfi - 2v/6) 2 - (5h* uSe - 2v/6) 2 = (h s vfi - Sh* uSe ) 2 + 
2 (Shl S9 — 2v/9) (h & ud — 5h* vSe ) and to bound from Lemmas 5.1-5.3, for example 
for the second term: 



j "w\K,e- 5 K,se\G a < c J w 1/2 g* uS9 \h 



8 

vfi 



< cO(5 



2w\ 



w 



l 2 G a = 0(5 



2u+U 
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ii) By a simple algebraic manipulation (and using that J go = 1) we can write 



(as a side note, this implies l}^ h Q < ^s,b A> s i nce > °f course, Jo(^o)co > 1 by 
a Cauchy-Schwarz type inequality). Thus observing that lim^oO) = +oo and 
lim^o Jo(h ) = were, in fact, intermediary results of the paragraph i) above, we 
obtain the claimed result for the ^-inefficiency. 



Proof of Theorem 3.1. Again let w(X) := [a s b ^(A)] 2 and let g := g 5 v6 . We easily 
see, from the expression of the equivalents of J w 1 ^ 2 and J w of Lemma 5.4, that 
it is sufficient to prove that J w(go/g) ~ {9 2u /6 2u ) J w. Let us first study this 
integral with the unaliased versions g* := g* u 5e and g^ := g* u 56o in place of g and 
g . Starting from the definition (1.1) of g*, we can write 

^o*(a)/(0 o V(a)) = (e 2 + (\/s)T +1/2 (o 2 o + (\/6) 2 r iu+1/2) = i + nm 

with 

F(A) := (l + (^-^)/(e 2 + A 2 )r +1/2 -l. 

Note that F is integrable over H. and has the same sign as Q 2 — 6\. By Lemma 5.3, if 
v > 1/2 then we obtain / w(X)F(X/S)dX = o(J w) and thus / w(6 2u g* Q / (6 2u g*)) ~ 
J w. The same result for the case v — 1/2 can be easily obtained (from the closed 
expression of w in this case). It remains to bound the effect of the aliasing: 

J w{g* Q /g*) - J w{g Q /g) < J w \{g* Q g - g Q g*) / g* g\ 

< c J w^Ug-gog*)/^ 

< c J w ^ 2 \g*-g \ + c J w l l 2 \(g-g*)gl/g*\ 

where the second inequality results from the obvious bound w 1 ^ 2 /bg < 1. Now the 
first term of this sum is clearly 0(5 2v ) from Lemma 5.1; this also holds for the 
second because of the boundedness of g%j g* which results from that of F(X/5). 
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